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^] ' This is the third paper of a series devoted to the study of the global properties 

\ of Joguet's sample of 79 nearby galaxies observable from the southern hemisphere, of 

• which 65 are Seyfert 2 galaxies. We use the population synthesis models of Paper II to 

derive 'pure' emission-line spectra for the Seyfert 2's in the sample, and thus explore 
l/^ , the statistical properties of the nuclear nebular components and their relation to the 

■ stellar populations. We find that the emission line clouds suffer substantially more 

extinction than the starlight, and confirm the correlations between stellar and nebular 

^ ^, velocity dispersions and between emission line luminosity and velocity dispersions, 

I ' although with substantial scatter. Nuclear luminosities correlate with stellar velocity 

O I dispersions, but Seyferts with conspicuous star-forming activity deviate systematically 

towards higher luminosities. Removing the contribution of young stars to the optical 
^ ' continuum produces a tighter and steeper relation, L oc erf , consistent with the Faber- 

Jackson law. 

Emission line ratios indicative of the gas excitation such as [OIIIJ/II/3 and 
[OIII]/[OII] are statistically smaller for Seyferts with significant star- formation, im- 
plying that ionization by massive stars is responsible for a substantial, and sometimes 

■ even a dominant, fraction of the Rf] and [Oil] fluxes. We use our models to constrain 

the maximum fraction of the ionizing power that can be generated by a hidden AGN. 
We correlate this fraction with classical indicators of AGN photo-ionization: X-ray lu- 
minosity and nebular excitation, but find no significant correlations. Thus, while there 
is a strong contribution of starbursts to the excitation of the nuclear nebular emission 
in low-luminosity Seyferts, the contribution of the hidden AGN remains elusive even 
in hard X-rays. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: starburst - galaxies: statis- 
tics 
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The emission-line fluxes and profiles of active galactic nu- 
clei (AGNs) carry, at least in principle, all the information 
necessary to study the properties of both the broad-line 
(BLR) and the narrow-line (NLR) regions, such as ioniza- 
tion and excitation, extinction, metallicity, kinematics, and 
AGN power (Heckman et al., 1981; Osterbrock & Shuder 
1982; Wilson & Nath 1990; Whittle 1992a,b,c; Nelson & 
Whittle 1996; Ho et al. 2003; Hao et al. 2005). 

An important and unavoidable issue in the study of 
the emission-line spectra, however, is contamination by the 
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underlying stellar populations, particularly for the Balmer 
lines. Historically, McCall, Rybski, & Shields (1985) recom- 
mended an empirical absorption correction of 2 A to the 
equivalent widths (EWs) for Ha, H/3 and H7, which, of 
course, is rather coarse. Nowadays, with available high spec- 
tral resolution stellar population synthesis models, such as 
Bruzual & Chariot (2003), or Gonzalez Delgado et al (2005) 
for example, it is possible to fit the observed absorption lines 
and emission-line-free continua simultaneously and thus to 
obtain the pure crnission-linc spectra after subtraction of the 
synthesized stellar components (Tremonti 2003; Kauffmann 
et al. 2003a,b; Hao et al. 2005; Cid Fernandes et al. 2005). 

This is the purpose of the present paper, the third of a 
series devoted to the study of the stellar populations in the 
nuclear regions of Seyfert 2 galaxies of Joguet's sample of 
nearby galaxies. The sample comprises 79 southern galaxies 
{6 < 15'^'*^) classified as Seyfert 2's with redshifts lower than 
z=0.017. A full description of the sample is presented in 
Joguet et al. (2001; Paper I). After a careful inspection of 
the data set including a revision of the data in the Uterature 
we found in Paper H that of the 79 sample galaxies, one is a 
type 1 Seyfert, 65 are bona-fide Seyfert 2's, 4 are LINER's, 
4 are starburst/HII rmclei, and the remaining 5 are normal 
galaxies showing no nuclear omission lines. 

We found in Cid Fernandes et al. (2004; Paper II) that 
the nuclear regions of the true Seyfert galaxies in the sam- 
ple present remarkably heterogeneous star formation histo- 
ries where young starbursts, intermediate age, and old stel- 
lar populations all appear in significant, but widely variable 
proportions. We also found that a significant fraction of the 
nuclei show a strong featureless continuum (FC) component, 
but that this component is not always an indication of a hid- 
den AGN; it could also betray the presence of young dusty 
starbursts (Storchi-Bergmann et al. 2000; Gonzalez Delgado 
et al. 2001 and Cid Fernandes et al. 2001). 

The present paper, the last of the series dealing with 
Joguet's (2001) sample, is devoted to the study of the emis- 
sion line properties of the 65 Seyfert 2 nuclei in the sample 
and their relation with the stellar components. In Section 2 
we measure the emission line fiuxes from the pure emission 
line spectra and derive the nebular extinction from Balmer 
decrements. Section 3 presents the relation between stellar 
and nebular velocity dispersion and Faber-Jackson relation 
for Seyfert 2 galaxies. The implications and discussion of the 
results are given in Section 4 and conclusions in Section 5. 



2 THE DATA 

By subtracting the synthetic stellar components we were 

able in Paper II to obtain remarkably clean pure emission 
line spectra. Examples of these 'residual' spectra were pre- 
sented in Figures 2-10 of Paper II. Using standard tools 
within IRAF ^ we measured the strengths of the most promi- 
nent emission features in the residual spectra, including line 
fiuxes and measurement errors, line widths, and equivalent 
widths using the continua from the synthetic spectra. The 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Extinction of the stellar component (Ay) as a func- 
tion of the nebular extinction (^ABalmer-^ derived from the Balmer 
decrements assuming Case B recombination and the standard 
reddening law. Only H/3 and H7 are used to compute A.y"''"^"'" . 
The dashed line represents x = y. The open symbols correspond 
to Seyfert 2s with direct evidence for hidden broad-line regions 
(HBLR) from spectropolarimetric observations. 

results are presented in Table 1 including estimates for the 
measurement errors derived from the photon statistics. Since 
many of our nuclei show faint, broad non-Gaussian wings, 
the widths were measured using two methods: fitting single 
Gaussiaus to the [OIII]A5007 line, and marmally measuring 
the widths at half-intensity. The two measurements agree 
to better than 10% so it is safe to interpret the line widths 
as true velocity dispersions. The velocity dispersion from 
the Gaussian fits and the associated errors obtained using 
{ngaussfit in IRAF) are given in the Table 1 where for com- 
pleteness we also tabulate the stellar velocity dispersions 
(cr*) from Paper II. While our population synthesis code 
(Paper II) does not provide estimates of the internal errors 
associated with the stellar velocity dispersions, a compari- 
son with published results from other groups using different 
techniques, indicates that the typical errors in our values 
are ~ 20km/ s. The absorption-corrected hard X-ray lumi- 
nosities (2-10 keV) collected from the literature and the ref- 
erences to the corresponding publications are given in the 
last two columns of Table 1. The lower-limits correspond 
to Compton-thick sources - column densities larger than 
10^*cm~^ - for which the direct hard X-ray continua are 
completely absorbed. 

2.1 Derived parameters 

Before deriving physical parameters from the emission line 
strengths we performed the standard sanity check of com- 
puting the ratio of [OIII]A5007 to [OIII]A4959 which, un- 
der all but the most extreme physical conditions, should be 
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Table 1. Emission line properties of Seyfert 2 galaxies 



Galaxy 


[Oil]" 

3727 


H5" 
4101 


H7" 

4340 


Hell" 

4686 


4861 


0[III]" 

5007 


<^[oni] 

km s~ 




kin 


erg s 


Ref 


ESO 103-G35 


135.3 ± 8.5 


10.1 ± 2.4 


18.2 ± 3.1 


4.8 ± 1.5 


44.8 ± 5.2 


430.3 ± 15.6 


259 


± 


8.1 


114 


42.79 


1 


ESQ 104-Gll 


45.7 ± 5.0 


6.9 ± 1.8 


12.3 ± 2.5 


7.6 ± 1.6 


26.7 ± 3.6 


279. S ± 12.4 


177 


± 


2.8 


130 






ESO 137-G34 


719.0 ± 19.2 


66.3 ± 5.7 


107.0 ± 7.4 


55.3 ± 5.1 


240.0 ± 11.1 


2541.0 ± 36.0 


239 


± 


3.8 


133 






ESO 138-GOl 


902.3 ± 21.5 


107.6 ± 7.4 


197.0 ± 10.1 


120.0 ± 7.9 


433.0 ± 16.1 


3693.0 ± 44.5 


136 


± 


1.8 


80 


>41.46'' 


2 


ESO 269-G12 


46.0 ± 4.9 


3.3 ± 1.3 


4.8 ± 1.6 




7.8 ± 1.9 


68.6 ± 6.6 


150± : 


11.3 


161 






ESO 323-G32 


80.1 ± 6.3 


10.7 ± 2.8 


15.3 ± 2.6 


13.7 ± 2.4 


36.8 ± 4.1 


460.6 ± 15.6 


132 


± 


1.6 


131 






ESO 362-G08 


48.0 ± 5.1 


15.6 ± 2.8 


22.5 ± 2.6 


13.3 ± 3.2 


38.4 ± 4.4 


286.2 ± 12.2 


161 


± 


2.3 


154 






ESO 373-G29 


186.6 ± 9.7 


40.0 ± 4.4 


65.0 ± 5.7 


25.1 ± 3.6 


142.3 ± 8.5 


751.0 ± 19.9 


78 ± 1.3 


92 






ESO 381-G08 


444.7 ± 15.2 


62.8 ± 5.9 


127.5 ± 8.5 


39.5 ± 4.8 


301.0 ± 13.1 


1942.0 ± 32.7 


206 


± 


3.8 


100 






ESO 383-G18 


80.0 ± 6.3 


9.4 ± 2.2 


18.2 ± 3.0 




42.5 ± 4.7 


360.9 ± 13.6 


105 


± 


1.6 


92 






ESO 428-G14 


910.6 ± 21.7 


93.1 ± 7.4 


187.4 ± 9.8 


126.8 ± 7.8 


397.0 ± 14.6 


6198.0 ± 52.3 


187 


± 


2.2 


120 


>40.64'' 


3 


ESO 434-G40 


118.8 ± 7.9 


28.9 db 3.6 


37.6 ± 4.1 


32.7 ± 3.8 


84.5 ± 6.8 


913.4 ± 22.0 


106 


± 


1.4 


146 


43.06 


1 


Fairall 334 


26.9 ± 3.7 


3.6 ± 1.3 


6.4 ± 1.7 


± 


14.3 ± 2.7 


86.6 ± 6.8 


130 


± 


2.8 


104 


42.95 


4 


Fairall 341 


40.0 ± 4.3 


7.7 ± 1.6 


9.7 ± 2.2 


8.8 ± 2.1 


25.1 ± 3.4 


279.2 ± 12.1 


120 


± 


1.6 


122 


43.03 


4 


IC 1657 


34.3 ± 4.2 




4.7 ± 1.5 


± 


9.7 ± 2.2 


25.7 ± 3.6 


148 


± 


6.5 


143 






IC 2560 


175.9 ± 9.6 


40.5 ± 4.5 


61.9 ± 5.6 


49.2 ± 5.0 


135.0 ± 8.0 


1250.0 ± 25.8 


135 


± 


1.6 


144 


>40.94'' 


5 


IC 5063 


247.4 it 11.4 


22.7 i 3.5 


41.5 i 4.7 


15.6 ± 2.6 


106.0 ± 7.5 


911.4 ± 22.2 


183 


± 


2.7 


182 


42.87 


1 


IC 5135 


305.2 ± 12.4 


34.6 ± 4.7 


80.0 ± 6.7 


32.8 ± 4.0 


191.6 ± 10.5 


1108.0 ± 24.6 


316 


± 


8.5 


143 


>41.40'' 


6 


F11215-2806 


51.8 ± 6.2 


4.6 ± 1.7 


8.1 ± 2.1 


6.6 ± 1.9 


22.5 ± 3.6 


243.2 ± 11.9 


184 


± 


2.4 


98 






M+01-27-20 


66.4 ± 6.7 


9.2 ± 2.1 


15.7 ± 2.8 


6.4 ± 1.7 


39.9 ± 4.5 


167.0 ± 9.0 


103 


± 


2.1 


94 






M-03-34-64 


310.3 ± 13.1 


62.8 ± 6.1 


130.4 ± 8.7 


146.3 ± 8.2 


320.9 ± 14.0 


4434.0 ± 47.2 


461 


± 


7.6 


166 


42.63 


1 


Mrk 897 


101.7 lb 7.1 


23.6 ± 3.7 


44.9 ± 4.7 




103.0 ± 7.1 


55.8 ± 5.6 


154 


± 


34.5 


133 






Mrk 1210 


212.1 ± 10.8 


55.1 ± 6.0 


131.9 ± 8.7 


64.3 ± 6.8 


282.0 ± 13.0 


2847.0 ± 40.9 


275 


± 


3.9 


114 


>41.96'' 


1 


Mrk 1370 


49.5 ± 5.1 


5.1 ± 1.5 


8.5 ± 1.9 




19.6 ± 3.2 


128.8 ± 8.4 


188 


± 


3.9 


86 






NGC 424 


455.2 ±16.7 


147.2 ± 9.1 


288.8 ± 12.9 


204.1 ± 10.8 


766.0 ± 21.1 


4267.0 ± 47.1 


298 


± 


5.3 


143 


>41.60'' 


2 


NGC 788 


99.2 ± 7.1 


10.4 ± 2.3 


17.3 ± 2.9 


8.8 ± 2.0 


31.6 ± 4.0 


364.8 ± 13.7 


101 


± 


1.5 


163 






NGC 1068 


3379.0 ± 43.0 


744.8 ± 20.5 


1106.0 ±28.3 


1102.0 ±25.9 


3194.0 ±43.3 


3.558e±04±135.4 


561 


± 


9.0 


144 


>40.98'' 


1 


NGC 1125 


86.3 ± 6.7 


7.9 ± 2.0 


16.0 ± 2.8 


6.3 ± 1.8 


35.1 ±4.3 


227.9 ± 11.3 


186 


± 


3.2 


106 






NGC 1667 


38.7 ± 4.6 




7.8 ± 1.8 




10.5 ±2.2 


91.3 ± 6.9 


174 


± 


3.1 


149 


>42.34'' 


1 


NGC 1672 


97.8 ± 7.2 


17.9 ± 3.4 


46.8 ± 4.8 




111.1 ±7.3 


83.3 ± 6.7 


129 


± 


42.3 


97 


41.46 


1 


NGC 2110 


944.5 ± 22.3 


55.6 ± 6.5 


98.5 ± 7.5 


19.2 ± 3.0 


167.1 ±9.5 


751.9 ± 20.7 


199 


± 


7.1 


242 


42.57 


1 


9Q7Q 


O'i.O ZIZ 


4 8 i 1 8 


7 5 ± 2 




11 6 ± 2 4 


1 no 7 -1- 7 


170 


± 


2.9 


112 






NGC 2992 


125.6 ± 8.1 


13.2 ± 2.8 


24.4 ± 3.9 


10.7 ± 2.2 


63.8 ± 6.2 


519.9 ± 17.4 


164 


± 


2.2 


172 


41.71 


1 


NGC 3035 


37.4 ± 4.3 




7.8 ± 1.9 




19.7 ± 3.6 


171.2 ± 9.5 


121 


± 


1.9 


161 






NGC 3081 


288.8 ± 12.2 


46.6 ± 4.8 


78.0 ± 6.4 


75.1 ± 6.0 


172.3 ± 9.6 


2109.0 ± 33.6 


93 ± 1.4 


134 


41.91 


1 


NGC 3281 


26.7 ± 3.6 


6.1 ± 1.6 


6.4 ± 1.8 


6.4 ± 1.6 


17.4 ± 2.9 


129.4 ± 8.2 


141 


± 


3.2 


160 


42.79 


1 


NGC 3362 


67.9 ± 6.6 


6.0 ± 1.7 


10.4 ± 2.6 


6.4 ± 1.9 


30.5 ± 3.9 


237.2 ± 11.5 


179 


± 


3.8 


104 






NGC 3393 


668.4 ± 17.3 


72.1 ± 6.0 


122.3 ± 8.1 


76.9 ± 6.1 


270.5 ± 12.0 


2680.0 ± 37.9 


184 


± 


2.4 


167 


>41.08'' 


1 


NGC 3660 


23.9 ± 3.6 


6.6 ± 1.6 


13.2 ± 2.9 


7.1 ± 1.6 


39.8 ± 5.2 


106.2 ± 7.4 


90 ± e 


i.3 


96 






NGC 4388 


601.9 ± 16.1 


30.6 ± 3.9 


56.4 ± 5.5 


31.9 ± 4.0 


145.1 ± 8.8 


1621.0 ± 29.6 


152 


± 


2.0 


111 


42.76 


1 


NGC 4507 


1022.0 ± 23.5 


134.6 ± 8.7 


262.0 ± 12.3 


100.4 ± 7.5 


578.8 ± 18.5 


4455.0 ± 51.7 


213 


± 


3.6 


144 


43.28 


1 


NGC 4903 


52.3 ± 4.8 


8.2 ± 2.2 


10.3 ± 2.4 


8.8 ± 1.8 


28.7 ± 3.7 


262.0 ± 11.6 


158 


± 


1.8 


200 






NGC 4939 


260.5 ± 11.2 


39.8 ± 4.4 


61.0 ± 5.6 


62.1 ± 4.9 


138.5 ± 8.6 


1699.0 ± 29.0 


163 


± 


1.8 


166 


41.96 


1 


NGC 4941 


276.5 ± 11.9 


27.7 ± 3.7 


47.5 ± 5.0 


36.7 ± 3.8 


114.4 ± 7.8 


1426.0 ± 27.9 


119 


± 


2.2 


98 


40.90 


1 


NGC 4968 


40.6 ± 4.7 


10.3 ± 2.3 


14.5 ± 2.7 


9.1 ± 2.3 


36.7 ± 4.5 


434.4 ± 15.5 


176 


± 


1.9 


121 


>40.87'' 


1 



equal to 3.0 (e.g. Rosa 1985). Using only the 62 objects for 
which [OIII]4959 is detected with good S/N we find a mean 
value of 3.00 ± 0.08. Therefore, we conclude that our data 
are not affected by unpleasant instrumental effects (such as 
saturation or non-linearities) which could compromise the 
reliability of the line ratios. 



In order to analyze the emission line spectra, and for 
lack of a better option, wo will assume that the nuclear neb- 
ulae of type-2 Seyferts are similar to giant HII regions and 
use the classical approach (e.g. Peimbert & Torres- Peimbert, 
1981). 



2.1.1 Extinction 

Assuming Case B recombination theory and a standard red- 
dening law (Cardelli, Clayton & Mathis 1989), one obtains 
the logarithmic reddening correction at H/3, C(H/3), from 

the Balmer decrements (in our data from H7 and H5; see 
e.g. Torres- Peimbert, Peimbert & Fierro 1989), which are 



C{H(3) = -7.41 X log( 



C{HI3) = -4.95 X log( 



iHf/lHg 

Fhs/Fhi3 ^ 

Ihs/Ihi3 



for 



for H6 



(1) 



(2) 



Where Fh-//Fh0, Fh5/Fhi3 and Ihj/Ihi}, Ihs/Ihi} are the 
observed and intrinsic BaJmer line ratios, respectively. In 
this paper we adopt the intrinsic ratios to be 0.466 and 



4 Gu, Melnick, Cid Fernandes, Kunth, Terlevich, Terlevich 



Table 1. Emission line properties of Seyfert 2 galaxies (continued) 



Galaxy 


[Oil]" 






HcII'^ 


H/3" 


OIIII]" 


°"[OIII] 




Lx 


R.cf 






3727 


4101 


4340 


4686 


4861 


5007 


km s~ ^ 


km s~ 


erg s^"^ 




NGC 


5135 


133.1 it 8.3 


21.0 i 3.2 


41.0 it 4.5 


20.6 it 2.8 


99.6 ± 7.3 


363.6 it 14.5 


155 it 4.4 


143 


>40.86 


1 


NGC 


5252 


186.6 ± 9.8 


12.2 ± 2.5 


21.7 it 3.4 


4.6 it 1.2 


40.8 ± 4.6 


268.1 ± 12.0 


195 it 6.9 


209 


43.07 


1 


NGC 


5427 


38.0 ± 4.5 


7.0 ± 1.9 


8.3 ± 1.9 


3.9 it 1.3 


20.0 ± 3.3 


205.4 ± 11.1 


161 it 2.1 


100 






NGC 


5506 


513.6 ± 16.2 


38.0 ± 4.7 


79.5 it 6.5 


28.3 it 3.8 


206.8 ± 10.7 


1614.0 ± 29.8 


203 it 3.0 


98 


42.89 


1 


NGC 


5643 


753.7 ± 19.8 


57.4 ± 5.6 


104.9 it 7.4 


61.7 it 5.7 


224.3 ± 10.8 


2393.0 ± 36.0 


154 it 2.4 


93 


>40.60'' 


1 


NGC 


5674 


21.6 ± 3.3 


3.3 ± 1.3 


4.7 ± 1.6 


1.2 it 0.7 


9.3 it 2.3 


62.7 it 5.9 


145 it 4.4 


129 


43.18 


1 


NGC 


5728 


253.1 ± 11.3 


29.0 ± 3.8 


50.0 it 4.9 


28.5 it 3.8 


124.9 ± 7.9 


1154.0 ± 25.0 


210 it 4.1 


155 






NGC 


5953 


121.4 ± 7.8 


17.6 ± 3.5 


37.7 it 4.7 


8.1 it 1.2 


74.8 ± 6.1 


222.2 it 10.9 


180 it 13.7 


93 






NGC 


6221 


201.5 ± 10.1 


83.6 ± 6.9 


168.1 it 9.4 


7.6 it 2.5 


385.7 ± 14.6 


249.5 it 11.7 


244 it 22.0 


111 






NGC 


6300 


55.3 ± 5.3 


4.9 ± 1.6 


6.0 ± 1.8 


6.0 it 1.8 


15.5 ± 2.8 


234.6 it 11.3 


121 it 4.7 


100 


41.23 


7 


NGC 


6890 


64.2 ± 5.9 


14.5 ± 2.5 


19.2 it 3.3 


19.1 it 3.2 


50.2 ± 5.2 


720.7 it 20.1 


204 it 3.0 


109 






NGC 


7172 


4.8 ± 1.6 








3.2 it 1.1 


15.5 it 2.8 


132 it 19.9 


190 






NGC 


7212 


476.3 ± 15.8 


49.5 ± 5.1 


100.8 it 7.3 


53.3 it 5.3 


242.1 ± 11.6 


2890.0 ± 39.1 


205 it 2.5 


168 


42.24 


6 


NGC 


7314 


22.0 ± 3.3 


3.7 ± 1.5 


8.5 ± 2.1 


6.8 it 1.8 


19.1 ± 3.3 


216.4 it 10.5 


100 it 1.4 


60 


42.22 


1 


NGC 


7496 


135.1 ± 8.3 


32.0 ± 3.9 


58.4 it 5.4 




150.3 ± 8.9 


96.4 it 7.3 


205 it 21.0 


101 


41.66 


4 


NGC 


7582 


140.7 ± 8.5 


49.9 ± 5.3 


98.3 it 7.2 


29.8 it 3.5 


276.0 ± 12.0 


570.1 it 17.2 


127 it 5.0 


132 


42.16 


1 


NGC 


7590 


55.1 ± 5.3 


7.6 ± 2.4 


9.8 ± 2.2 


5.9 it 1.4 


24.0 ± 3.3 


112.2 it 7.5 


73 it 8.7 


99 


40.81 


1 


NGC 


7679 


147.0 ± 8.7 


24.0 ± 3.4 


48.0 it 4.8 


8.0 it 2.8 


114.0 ± 7.8 


154.5 it 9.9 


176 it 11.7 


96 


42.47 


8 


NGC 


7682 


262.0 ± 11.6 


18.4 ± 3.0 


32.6 it 4.2 


17.8 it 2.9 


71.6 ± 6.3 


729.2 it 20.0 


139 it 2.3 


152 


42.86 


4 


NGC 


7743 


83.1 ± 6.7 


7.6 ± 2.2 


8.5 ± 2.3 




15.7 ± 2.8 


79.2 it 6.6 


160 it 31.0 


95 


39.51 


9 



" emission line flux in units of 10 erg cm ^. Compton-thick source with column density > lO^^cm ^. Notes: sources of X-ray 
luminosity - (1) Bassani et al. 1999; (2) CoUinge & Brandt 2000; (3) Maiolino et al. 1998; (4) PoUetta et al. 1996; (5) Iwasawa, Maloney 
& Fabian 2002; (6) Risaliti et al. 2000; (7) Matsumoto et al. 2004; (8) Delia Ceca et al., 2001; (9) Terasliima et al. 2002. 



0.256, respectively (Osterbrock 1989). One then compares 
the two values which for properly exposed and calibrated 
data should not differ by more than a few percent. In our 
case, even after subtraction of the stellar component, H5 
is detected with good S/N for 42 objects, out of which 
9 have unphysical Balmer ratios (i.e. Fhs/Fhp > 0.256). 
H7 are detected with good S/N for all objects except NGC 
7172, and 14 Seyfert 2s have unphysical Balmer ratios (i.e. 
Fhi/Fhp > 0.466). For 33 objects with both good S/N 
and physical Balmer ratios for H7 and YLS, we obtain an 
average difference between the two Balmer decrements of 
AC(H/3)= 0.02 ± 0.27(a), which indicates that it is prob- 
ably safe to use the Balmer decrements to infer the value 
of the extinction for these objects. In the following of this 
paper, we will simply use Fh-i/Fhi3 and Equation (1) to cal- 
culate C(H/3) for Seyfert 2 galaxies and exclude NGC 7172 
from further analysis. For objects with unphysical Balmer 
decrements we will assume C(H/3)=0. 

We can investigate this further by comparing the nebu- 
lar extinction to the stellar extinction (A^) from our popu- 
lation synthesis models. FigureQshows the run of Ay versus 
j^Baim<,r ^ g gg^ ^ 2.5 X C (H (3) . The nebular extinction is 
always larger as expected since the stellar population is a 
mix of unredenned old stars with reddened young stars in 
varying proportions (Calzetti, Kinney & Storchi-Bergmann 
1994; Mas-Hesse & Kunth 1999). The open symbols repre- 
sent objects where the spectropolarimetric observations de- 
tected the polarized broad emission lines, mainly from the 
compilation by Gu & Huang (2002) plus new observations 
by Lumsden, Alexander & Hough (2004; see also Paper II). 
In what follows we will correct our emission line data for 
extinction using the Balmer decrements, but we note that 
none of our conclusions actually depend on the reddening 
corrections. 




l°g <^[oiii] (km s 1) 

Figure 2. The stellar velocity dispersion cr« derived from the 
population synthesis models as a function of the velocity width 
of the emission lines obtained from the [OIIIJ5007 line profile. The 
uncertainty in (t« is typically 20 km s~^, and the measurement 
error in ct^qjjj] is estimated by a single Gaussian fitting. The 
solid line shows a ordinary least square (OLS) bisector fit to the 
data of slope 0.84±0.08 and the Spearman rank-order correlation 
coefficient (Rs) is 0.26. Symbols have the same meaning as in Fig. 
1. 
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Figure 3. Luminosity - velocity dispersion (L - a) relation for 
the nuclear regions of Seyfert 2 galaxies. The correlations of 
[OIIIJA5007 luminosity with the width of the nebular emission 
lines (ctjoiii] - left panel) and with the stellar velocity dispersion 
CT* (right panel) are shown. Symbols have the same meaning as 
in Fig. 1, and the lines show ordinary least square (OLS) bisector 
fits to the data. 



Figure 4. The relationship between monochromatic continuum 
luminosity {L'^^^) and stellar velocity dispersion for Seyfert 2 
galaxies. Upper panel: the continuum luminosity is estimated 
from the total 4861 A continuum emission; Lower panel: Con- 
tinuum luminosity is estimated using only the luminosity from 
the old contribution. Squares (triangles) indicate Seyfert 2s with 
Wjf < 10A( (> 10 A). The y-axis are in units of ergs~^A~^. 



3 RESULTS 
3.1 Kinematics 

Our stellar population synthesis code uses a Gaussian broad- 
ening function to match the synthetic stellar populations to 
the observed galaxy spectra, and thus provides reliable es- 
timates of the stellar velocity dispersion cr* for each galaxy. 
We showed in Paper II that our dispersions are in good 
agreement with results obtained using other methods and 
by other authors (mostly Nelson & Whittle 1995; but see 
also Garcia-Rissmann et al., 2005). We have used the re- 
sults of Paper II to correct the emission line widths for 
instrumental broadening as cr[oiii] = v^o()s^^^ist' ^here 
finst = 62 kms~^. In what follows we will always use the 
corrected velocity dispersions in our analysis of the data. 

FigureElpresents a plot of (corrected) stellar versus neb- 
ular velocity dispersions. The solid line is the best fit to the 
data by means of the ordinary least squares (OLS) bisector 
method (Isobe et al. 1990), logcr. = (0.84±0.08) logcr[oiiii + 
(0.24 ±0.17). Figure |5|also quotes the Spearman rank-order 
correlation coefficient {Rs, Press et al. 1992) to be 0.26, and 
a probability of Pnuii = 0.03 for the null hypothesis of no 
correlation between a* and (T[oiii]. Our data thus confirms 
the findings of Nelson & Whittle (1996) and Jimenez-Benito 
et al. (2000) of a correlation between nebular and stellar 
velocity dispersions, albeit with a very significant scatter. 
In general we expect the emission line-profile widths to be 
due to the combination of hydrodynamical effects (winds, 
jets) and the motions of individual line emitting clouds in 
the gravitational potential of the nuclear regions. Thus, if 



the nebular and the stellar components feel the same grav- 
itational field, then the emission lines should be systemati- 
cally broader than the stellar lines. In fact, 18 of our Seyfert 
2's show stellar velocity dispersions larger than the nebular 
ones, although only for 7 of these, (ESO 434-G40, NGC 788, 
NGG 2110, NGG 3035, NGG 4903, NGC 7172 and NGC 
7590), the differences are larger than the typical errors in 
stellar velocity dispersion (20km/s). For NGC 7172, the dif- 
ference is due to the low signal-to-noise of the spectrum, 
while 4 of the remaining 6 have clearly non- Gaussian emis- 
sion lines with real widths significantly broader than the 
single Gaussian fits. Thus, for only two objects in our sam- 
ple - NGC 788 and NGC 3035 - the stellar lines appear to 
be significantly broader than the nebular lines. 

It is well known that [OIII] luminosity is a good tracer 
of the strength of the nuclear activity (Mulchaey et al. 
1994; Kauffmann et al. 2003b), Figurej^shows the relations 
between I/[oiii] and the velocity dispersions of the nebu- 
lar and stellar components. We remark that, although our 
extinction-corrected [OIII] luminosities have been observed 
and reduced in a homogeneous way, there is still an intrinsic 
source of scatter arising from the fact that, as a function of 
distance, our spectrograph slit encompasses different frac- 
tions of the total nuclear emission for different objects. The 
lines show ordinary least square (OLS) bisector fits to the 
data of slopes 3.8 ± 0.4 if we use o"[oiii], and 4.1 ± 1.0 for 
a,. The Spearman rank-order correlation coefficients, Rs, 
are 0.46 and 0.34 and Pnuii < 10~* and 0.007, respectively. 
Two features of these plots will be important in the dis- 
cussion: (1) the correlation is clearly better for the nebular 



6 Gu, Melnick, Cid Fernandes, Kunth, Terlevich, Terlevich 




5 10 15 20 



Wk (A) 

Figure 5. Relation between gaseous excitation (Upper: 
[OIII]/H/3 and Lower: [OII]/[OIII]) and equivalent width of Call 
K (a good indicator of the age of nuclear stellar populations) from 
population synthesis models. Symbols have the same meaning as 
in Fig. 1. Objects with HBLRs (open symbols) are not systemat- 
ically different from the rest of the sample in these plots. 



velocity dispersions, and (2) the objects with indications of 
hidden broad line regions lie predominantly above the best 
fit line in the L[oiii] — cr[oni] plot. 

The Faber-Jackson relation (Faber & Jackson 1976) 
is the correlation between bulge (continuum) luminosity 
{LbuIqe) and central (stellar) velocity dispersion. Since we 
do not have L Bulge, we will use the monochromatic contin- 
uum luminosity at A = 486lA, LJf^*, as a proxy. Figure 2] 
(upper panel) shows the relation between the I/^^' and cr*. 
We have used different symbols to parametrize the nuclear 
stellar populations by means of Wk - the equivalent width 
of the Call K absorption line - which is a very good indicator 
of the mean age of the stellar population (e.g. Cid Fernandes 
et al. 2001; and Paper II). Thus, triangles represent young 
objects ( Wk < lOA) and squares galaxies with Wk > lOA 
corresponding to objects whose continuum is dominated by 
old populations. The plot is quite scattered, but shows inter- 
estingly (top panel) that Seyfert 2s with Wk < lOA lie sys- 
tematically above the best fit line, indicating lower M/L ra- 
tios. If we use the continuum from the old populations only, 
L^^' X Xoid, where Xoid is the fraction of stars older than 
10^ yrs from our stellar population synthesis models (Pa- 
per II), the scatter is significantly reduced and the 'young' 
objects are more symmetrically distributed around the best 
fit line (lower panel). Furthermore, the slope changes from 
3.0 ± 0.5 (top panel) to 4.7 ± 0.5 (lower panel), closer to 
the actual Faber-Jackson relation. We confirm, therefore, 
that starburst activity distorts the Faber-Jackson relation of 
Seyfert galaxies as originally proposed by Nelson & Whittle 
(1996) and Gu et al. (1998). 



3.2 Physical conditions of the nebular gas 

The excitation of the nebular gas provides us with a fur- 
ther test of the link between the source of ionizing pho- 
tons and the nuclear stellar populations. Figure |3 plots 
two indicators of nebular excitation ([OIII] A5007/H/3 and 
[OII]A3727/[OIII]A5007) as a function of Wk described 
above. The Spearman tests gives Rs 0.22 and -0.58 (Pnuii 
0.09 and < 10^*), respectively. The trend of younger nu- 
clei having lower excitation and vice versa is clearly seen in 
[Oil] /[OIII], but the separation is less clear in [0III]/H/3. 
This is probably due to circum-nuclear starburst activity 
that powers a substantial fraction of the II/3 emission, but 
little of the [OIII] and Hell fiuxes in nuclei with low Wk- 
This would result in a dilution of the line ratios relative 
to the more extreme values attained by AGN without con- 
spicuous star-formation (large Wk', c.f. Cid Fernandes et al. 
2001). The fact that [OII]/[OIII] shows a clear trend suggests 
that [Oil] has a significant contribution from circumstellar 
HII regions. A similar result was found in composite low lu- 
minosity QSO spectra from the 2dF/6dF surveys (Croom et 
al 2002; but see Ho 2005). 

There is also a relation between the nuclear population 
ages and H/3 luminosity, as shown in Figure [S] which also 
shows that younger objects are more luminous in H/3. It can 
be verified using the data in Table 1 that these trends hold 
even at fixed velocity dispersions. That is, for a given value of 
(T* , younger objects tend to be systematically above the cor- 
relation line shown in Figure |3 thus indicating that trends 
in luminosity with age refiect real changes in the stellar mix. 



4 DISCUSSION 

4.1 Starbursts as sources of ionizing photons 

The young stellar components detected in many Seyfert 2s 
are bound to contribute to the ionization of the surround- 
ing gas. One should thus expect a proportionality between 
the optical continua associated with these components and 
the H/3 luminosity. This is confirmed in Figure [7| where we 
plot the observed H/3 emission line fluxes as functions of the 
ionizing stellar and featureless continua at 486lA, both cor- 
rected for the extinction derived from the synthesis models 
of Paper II (see also Cid Fernandes etal. 2001). The lines 
correspond to constant equivalent widths of lOA, lOOA, and 
lOOOA in each plot. From the left to the right panels of Fig- 
ure 7, the Spearman ranks are 0.60, 0.71 and 0.84 respec- 
tively showing that the correlation gets much better when 
we only use the optical flux associated with the putative 
ionizing continua (xY/Fcfn/?') instead of the total stellar 
continuum light. 

We can also use our stellar population synthesis results 
to predict the H/3 luminosities from the number of Lyman 
continuum ionizing photons, Qh, emitted by each stellar 
component using the theoretical values given by Bruzual & 
Chariot (2003) for each stellar base. Since the fiux of Ly- 
man continuum photons, Qh, changes by more than one 
order of magnitude during the first few Myrs of a starburst, 
whereas the optical spectrum changes relatively little, the 
young (5 Myr) component of our models represents the opti- 
cal light very well, but the corresponding Qh may be wrong 
by factors of 10 or more. Therefore we have calculated an 
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Figure 6. Relation of H/3 luminosity as a function of 3 indicators of the mean age of the nuclear stellar population: Wk, mean age 
(log t), and fraction of the total continuum emitted by young stars to the featureless (non thermal) continuum Xy/FC- Symbols have 
the same meaning as in Fig. 1. 



extended grid of models adding two additional young com- 
ponents of 1 Myr and 3 Myr. 

FigurelHIa) compares the observed Lh/s with the values 
predicted by our models using a single age of 5 Myr for 
the young stellar component, LJ^^. No clear correlation is 
seen, and in fact the observed luminosities are systematically 
larger than the predicted ones (the one-to-one correlation 
is indicated by the dashed line). As discussed above, 
depends strongly on the ages of the youngest components 
included in the synthesis models. This is illustrated in Figure 
ISJc) which shows the comparison between the predicted and 
observed luminosities obtained using the extended models. 

The predicted luminosities, however, (Figure |HJ;) still 
do not include the FC component in the computation of 
Qh- Figure |7| demonstrates that, regardless of its true na- 
ture (young stars or reflected hght from a hidden AGN, 
Storchi-Bergmann et al. 2000; Gonzalez Delgado et al. 2001 
and Cid Fernandes et al. 2001), the FC component is as- 



sociated with ionizing photons, so it must be included in 
the calculation of Qh, the question is how. Here we have 
converted the FC continuum luminosities from the models 
[xfcL'^h]}*) to emission-line luminosities ( L^^) by assum- 
ing that photo-ionization by the FC component alone would 
result in an EW(H/3) of lOOA, typical value for both star- 
burst and Seyfert nuclei (Figure |7|l. This estimate is then 
added to L^p to compute the total predicted H/3 luminos- 
ity. The results are shown in panels (b) and (d) for the two 
stellar bases considered above. In both cases we obtain very 
good correlations, in particular when we use multiple young 
stellar components (d) for which the points are very well 
fitted by a line of slope one. 

In order to check the consistency of the models we have 
also plotted in Figure |H| the predictions for the 4 objects 
in Joguet's original sample that are pure starbursts. These 
objects, shown as stars, lie exactly where the predicted and 
observed H/3 luminosities are equal, thus validating the in- 
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Figure 7. Relation between the continuum at H/3 and H/3 emission line fluxes. The panels differ in which component of the continuum 
is used in the abscissa, (left) total observed continuum, (middle) only the continuum due to young components; (right) the continuum 
from young and FC components. All continua are given in units of ergs~^cm~^A~^. The dashed lines represent, from left to right, loci 
of constant equivalent widths of 10, 100 and 1000 A. The objects with HBLRs (open symbols) are not segregated in these plots. 



ternal consistency of our models. It is interesting to notice 
that, while for 3 of these objects the synthesis yields a rela- 
tively small FC components (due to the younger starbursts) , 
for one, NGC 3256, the inclusion of the FC increases the pre- 
dicted H/3 luminosity by a factor of ~ 3, which illustrates 
the ambiguity of the FC discussed above and in Paper II. 

We can now use our models to calculate the fraction of 
the ionization flux that is produced only by the starburst 
components, fsB = L'^'^/L^^'^, shown in Figure |n| where 
for i^^"^ we have summed all ionizing photons supplied by 
the stellar components in the base except for the FC compo- 
nent, though it could be due to a young dusty starburst. The 
median fraction of starburst ionizing power is 65%, while 
for only 20% of the sample the starburst contribution is less 
than 10% of the total power required to ionize the gas. These 
figures are lower-limits since with our data we cannot dis- 



tinguish AGN from dusty starbursts, but, even if all the FC 
came from the AGN, the distribution is still consistent with 
50% of the objects being fully powered by starbursts. 

To summarize, besides the presence of weak (scattered) 
broad-line components in some nuclei, our optical spectra do 
not provide us with reliable direct tracers of the properties, 
nor even of the existence, of the hidden AGN in the nuclear 
regions of nearby low-luminosity Seyfert 2 galaxies. 



4.2 Starbursts versus Monsters 

Within the ambiguities of the FC discussed above, the pa- 
rameter fAGN ~ i — fsB should provide us with at least 
a hint of the ionizing power of the elusive hidden AGN. 
Figure [TI71 plots the relation between fAGN and nebular ex- 
citation. In order not to crowd the plots we have left out 
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Figure 8. Predicted versus observed H/3 luminosities (all in ergs~^). L'^g = 4.76 X 10~^'^Q^ is the H/3 luminosity powered exclusively 
by the stellar components in the base. See text for the FC contribution to the rate of ionizing photons, (a) L'^^ versus L"^^ for the 5 
Myr stellar base used in this paper, (b) Same as (a), but adding to L*jjp- (c) obtained with a larger stellar base (1, 3 and 5 

Myr). (d) As (b), but for the enlarged base used in (c). As usual HBLR objects are plotted with open-squares. The 4 starburst galaxies 
in Joguet's sample are shown as open stars. 



the Hell lines for which we do not see any significant cor- 
relations, and which may be compromised by Wolf-Rayet 
features. 

While these plots show an intriguing structure, there is 
no obvious correlation between AGN fraction and excitation. 
Moreover, objects with putative HBLRs (open symbols) do 
not seem to stand out in these plots either. Since excitation 
is one of the classification criteria for Seyferts, these plots 
indicate that Jagn is not a very good starburst/monster 
discriminator. Another potentially powerful indicator of ac- 
cretion power is hard X-ray luminosity. Figure 1111 (left 
panel) shows the relation between fAON and absorption- 
corrected hard X-ray luminosity. The lower-limits (arrows) 
correspond to Compton self-absorbed objects (hydrogen col- 
umn densities > lO^'^cm"^) for which only refiected X-rays 



are observed. The right panel shows Lx as a function of 
[OIII]/[OII]. 

There is no correlation between hard Ljc and AGN frac- 
tion, but there is a clear trend of nebular excitation being 
higher for more X-ray luminous objects. These plots also 
show very clearly that nuclei with HBLRs have high X-ray 
luminosities, although the interpretation of this result is not 
straightforward since there may be a substantial contribu- 
tion to the X-rays from off-nuclear sources (Jimenez-Bailon 
et al., 2003; Levenson et al., 2005). It is therefore relevant 
to examine the relation between X-ray luminosity and the 
global parameters of the sample: velocity dispersion and H/3 
luminosity. These relations are presented in Figure lT^ where 
we see a weak trend of Lx with cr[oni] , but a significant cor- 
relation with L(H/3) which is not a distance effect because it 
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Figure 9. Histogram distribution of the ratio of the observed 
I^Hi3 to the values predicted by the stellar components in the 
base, (a): only one at 5 Myr; (b): 1, 3 and 5 Myr models. 



Figure 11. Relationships between fAGN, nebular excitation 
and hard X-ray luminosity. Arrows show Compton self-absorbed 
galaxies for which the X-ray luminosities are lower-limits. Galax- 
ies with direct indications of a hidden BLRs are plotted as open 
squares. 
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Figure 10. Relationships between the maximum Lyman ioniza- 
tion power produced by a hidden AGN component, /agjVi and 
two indicators of the excitation of the nebular gas. As before, 
open squares correspond to objects with HBLRs. 



is also present in the fluxes. This is shown in Figure lTsl where 
for clarity we have omitted the lower- limits. The line shows 
the OLS fit to the data of slope 1.1 ± 0.2 {Rs = 0.19) so, 
with large scatter, the X-ray power is roughly proportional 
to the Lyman ionizing power. Since we have no control over 
the spatial regions covered by the two data sets, we should 
be very careful about how to interpret this plot. Taken at 
face value it means that the bulk of the X-rays are produced 
by the same sources that photoionize the gas. This is consis- 
tent with the fact that the few objects with [OIII]/H/3 < 3 
(cf. Figure 10) lie below the best fitting line. 

Since even using deep high resolution X-ray (Chandra) 
and UV (HST) images it remains extremely difficult to see 
the AGN in nearby Seyfert 2's ( e.g.. Jimenez-Bailon et al. 
2003; Levenson et al. 2005), and since even when it is seen, 
there is still a significant contribution from circum-nuclear 
sources, the only way to answer the Starburst or Monsters 
question is to have deep high spatial resolution images of a 
substantial sample of nearby Seyfert 2 galaxies over a wide 
wavelength range. Until this is available, the interpretation 
of integrated observations will be plagued with the same 
uncertainties that we have confronted in the present work. 

The INTEGRAL satellite has recently unveiled a puta- 
tive new class of powerful, highly obscured. X-ray sources 
associated with massive stars: High Mass X-ray Binaries 
(HMXBs; Chaty and Filhatre, 2005). If HMXB's represent 
a stage in the evolution of high-mass binaries, they could be 
very abundant in the circumnuclear starbursts presented in 
many Seyfert 2s, and thus potentially explain their X-ray 
morphologies. 
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Figure 12. Correlations between (2-lOkeV) X-ray luminosity, 
nebular velocity dispersion (<J[oiii]), and H/3 luminosity. As be- 
fore, arrows correspond to Compton self-absorbed lower limits, 
and open symbols to galaxies with hidden BLRs. 




Figure 13. 2-lOkeV X-ray flux against H,8 flux. We have omitted 
the Compton self-absorbed sources for clarity. The line shows the 
best fitting line of slope unity. 



ionized gas, and confirm the existence of correlation between 
stellar and nebular velocity dispersion. 

(iii) Seyfert 2 nuclei follow a Faber-Jackson like cor- 
relation between nuclear continuum luminosity and stel- 
lar velocity dispersion, but nuclei which also harbor star- 
formation deviate systematically from this relation. Remov- 
ing the young stars from the continuum luminosity (using 
the synthesis results) strengthens the correlation substan- 
tially, and results in a slope more compatible with the actual 
Faber-Jackson relation. 

(iv) There is a relation between ionized gas velocity 
dispersion and emission-line luminosities. We also find that 
Seyfert 2s with indication of broad lines cither from the pure 
emission line spectrum or spectropolarimctric observations 
show larger [OIII] luminosities. 

(v) We confirm our earlier inference that photoioniza- 
tion by young stars can explain a substantial fraction of the 
nuclear emission-line luminosities of Seyfert 2 nuclei, but 
neither the parameters derived from our optical spectra, nor 
the hard X-ray luminosities provide accurate indicators of 
the AGN contribution. 

(vi) Albeit with substantial scatter, there is a correla- 
tion between hard (2-20kcV) X-ray luminosity and H/9 lu- 
minosity for our sample of low-luminosity AGN. The recent 
discovery of a new population of highly obscured luminous 
X-ray sources associated with high mass stars may yield 
some new clues to the interpretation of the properties of low- 
luminosity Seyfert galaxies. The off-nuclear X-ray sources 
detected by Chandra in some of these galEixies may, in the 
end, be powered by the same starbursts that are producing 
the bulk of the optical emission. 
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5 CONCLUSIONS 

We have presented a study of pure emission line spectra of a 
sample of 65 Seyfert 2 galaxies obtained by subtracting the 
synthetic stellar component from the observed spectra. Our 
main results are summarized as follows. 

(i) For most Seyfert 2 galaxies, the gaseous nebular ex- 
tinction deduced from the Balmer decrements appears to be 
systematically larger than the stellar extinction. 

(ii) We compare the velocity dispersion of the stars and 
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